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I. Introduction 

It is well known that an electron in an isolated Bloch band will oscillate with frequency 
eFa/h in the Brillouin zone when it is subject to a uniform and constant electric field F, 
where a is the lattice constant. It was first pointed out by Wannier that such oscillatory 
motion is associated with a ladder like spectrum, now commonly known as Wannier-Stark 
ladders. Wannier further proposed that sequences of such ladders exist even when the 
interband coupling is taken into account. [1] 

This proposal was immediately challenged by Zak, who pointed out that sharp Wannier-| 
Stark ladders cannot exist in a system with any reasonably continuous periodic potential, 
because the spectrum is absolutely continuous. [2] A numerical study on the density of 
states p{E) for a system with a cosine potential in an electric field also showed no evidence 
of the ladder structure [3]. Since the existence of sharp Wannier-Stark ladders in the one 
band approximation was rigorously established [4], it was then suggested that a correct 
account of the interband transitions will smear out the discrete spectrum [5]. Further, it 
has also been proved that an isolated group of n Bloch bands can give rise to n sequences 
of sharp ladders, so one must consider an infinite number of coupled bands in order to 
obtain a continuous spectrum. Recently, by decomposing the linear electric potential into 
a sawtooth part plus a staircase part, Emin and Hart claimed that the bands may be effec- 
tively decoupled because the interband matrix elements of the staircase potential is zero. 
This would give a proof for the existence of a sharp ladder spectrum, but the arguments 
leading to the vanishing of interband matrix elements have been critcized[6]. 

On the other hand, experimental and theoretical evidence in support of the existence 
of Wannier-Stark ladders also increased during the same period of time [7]. Krieger and 
lafrate clarified the controversies concerning the use of periodic boundary condition and the 
neglect of interband transitions. [8] Although they did not discuss the energy spectrum at 
all, they can show that the transitions of electrons induced by an incident electromagnetic 
field are allowed only if /S.E = eFa, as if these levels exist. In recent years, thanks to 
technological advancement on superlattices, the existence of a ladder structure seems to 
be firmly established. [9] 
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A clear picture about the nature of the W-S ladders should have emerged from the 
many years of contraversial debates: the ladders are not infinitely sharp levels, but they are 
resonant levels whose life times can be long under conditions of small intcrband couplings. 
In earlier numerical works, such resonant levels are located by searching for poles of the 
S matrix or the resolvant in the complex energy plane, with the level widths given by 
the imaginary part of the poles. [7] In this paper, we reveal and display the W-S ladder 
structure in the local density of states p{x,E), which can be calculated in a simple and 
direct way. As we will show in the main text, the W-S ladders appear naturally and very 
clearly in a plot of p{x, E) as sequences of mountains, whose widths in the E direction 
represent the level widths, and whose widths in the x direction represent the localization 
of the ladder states. Such structures can be easily washed out in the total density of states 
p{E). 

In order to simplify the calculation, we use in our model a chain of delta function 
barriers to simulate the periodic potential, and use a staircase potential to simulate the 
external field. There has been criticism concerning the use of the staircase potential, 
arguing that it will lead to vanishing interband transitions [10], however, such criticism is 
invalid as has been shown by Kleinman[6] and others. Although our model is used mainly 
for qualitative studies, it can also be used to simulate real solids when the spatial extent 
of the wave functions are large compared to the lattice constant as has been demonstarted 
by Bentosela et al in a similar situation. [7] 

The organization of this paper is as follows: In Sec. 2 we describe our model and derive 
a formula for the local density of states. In Sec. 3 we present the numerical results for the 
local and total density of states with a variety of parameter values. In Sec. 4, the behavior 
of the level width as a function of system size and other parameters is explained in a 
picture of wave tunneling through the potential barriers. 



II. The Model 
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We will study the wave function solution of a one dimensional model with Hamiltonian 



^2 ^2 

S <^(^-ja) + Kxt (1) 

where V^xt is a staircase potential, simulating an external electric field. Reference of energy- 
is chosen such that 



Vext = -jA for ja<x< {j + l)a , j e Z (2) 

We will take /2ma? as the unit of energy, therefore A = 0.2 corresponds to an electric 
field of F = 7 X 10^ (V/cm) for o = 200A and m = 0.067me. 

We will choose our boundary such that the potential beyond both ends of our crystal 
levels off, as shown in Fig.l. The cells are numbered by the integer j in (2). The wave 
function in cell j is 

iPjix) = aje'''^'= + bje-'''^'' (3) 

where kj = {7r/a)y/E + jA (new unit). The coefficients of cell j + 1 and cell j are related 
by a transfer matrix Mj : 



where 

_ 1 / {kj + kj+i - ifi)e-'ikj + i-k,)(j + l)a ^ /.^.^^ _ + X 

^ ~ 2kj+i V {-kj + kj+i + i«)e^(fc.+i+fc.)0+i)« (^kj + kj+i + iK)e^(fc.+i-fc.)0+i)« J 

(5) 

and K = TT^Vo/a. The above formula remains valid when E < —jA, in which case we take 
kj = —i{^^/a)^y\jA + E\. Notice that detM^ = kj/kj+i, which is not equal to the usual 
value of unity, because the velocities corresponding to a given energy in neighboring cells 
are different. 

In our study we will take E < NiA, and drop off the exponentially large solution at the 
left edge by letting 6_7V;-i = 0. The two solutions a^v^e*'^^'-^, bN^e~^^^^^ beyond the right 
edge correspond to incoming and outgoing waves. The normalization of wavefunctions is 
set by taking |ajv^| = 1. Furthermore, in the actual calculation, the number of cells on the 



negative x axis is limited to 10 for convenience. The error involved is negligible as long as 
we are considering energy only a few A's large. 

The quantity we are after is the local density of states defined by 



where Ek = {ka/ir)^ Unlike the total density of states p{E) = J dxp{x, E) which is usually 
calculated, the local density of states also contains information about the location of states, 
and is more appropriate for the purpose of revealing the nature of the W-S ladders. 



The coefficients Oj, hj[j — —{Ni + 1) ~ Nr) are calculated starting from 6_jv,-i = 0, 
and a_jv,-i is chosen such that lajv^j = 1. We can get p{x,E) for some fixed value of 
energy E. The energy is then changed to a different value and the calculation is done all 
over again. The result is double checked by calculating a-, , 6-, 's starting from oat^ and bN^. 

When the electric field is zero, we should see an energy band structure, as shown in 
Fig. 2. There are 50 cells (A^; = 10, A^^ = 40) in this system. Only the first band and 
the lower half of the second band are plotted. The band edges coincides very well with a 
simple theoretical calculation. 

Fig.3 is a plot of p{x, E) after we turn on the electric field. The lattice is composed 
of 50 cells ( labelled —10 ~ 39). Only 12 cells (—1 ~ 10) are plotted here. The energy 
range covered is 3 A's wide. Two sequences of states can be seen. The first sequence 
lies at energies equal to (0.34 + n) A, where n =integer and is equal to 0,1 and 2 in the 
energy range of the figure. The ladder at £■ = 0.34A has a primary peak in cell 1 and 2, 
with smaller peaks in the cells to the right. The second sequence is weaker and broader, 
but is still discernible, and lies at energies equal to (0.65 + ?i)A with n = 0, 1,2 in the 
figure. The ladder at -E = 0.65A is primarily peaked in cell 8, with the heights of the 
peaks slowly decreasing to the left and right. It is evident that p{x, E) is invariant under 
changes a;— >a; + a,£^— — Aasit should be. Each state in both sequences spans many 




(6) 



III. Result 
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lattice sites, and will not 'feel' much difference whether the electric potential is linear or 
staircase- like. 

The total density of states p{E) is obtained by integrating p(x, E) with respect to x 
up to the right edge of the lattice. The result is shown in Fig. 4 (a) (note the difference 
from the result in the absence of the field shown in Fig. 2). The two peaks within each 
A in energy correspond to the two sequences that can be observed in the local density of 
states. As we enlarge the system, the number of peaks increases (see Fig. 4 (b),(c),(d) ), 
corresponding to the increased number of observable sequences. One can imagine that as 
the system becomes infinite, there will be infinite number of peaks covering up the whole 
energy axis, making the ladder structure unrecognizable in p{E). Fig. 5 is a plot for an 
electric field that is five times larger than in Fig.4, with more peaks showing up in p{E) 
corresponding to more sequences of observable ladders in the system. 

For a system with {V,A,Nr) — (1,0.2,40), two sequences are observed, as in Fig. 3, 
but with much narrower level widths. In Fig. 6, the spatial extent of the states in these 
two sequences are plotted. One is at £^ ~ 0.3 A, the other is at ~ 0.7A. The state in 
Fig. 6(b) is essentially confined within 3 ~ 5 cells, which is more localized than the states in 
the first 'band' in Fig. 3. A cross section of it at x = 3.5a is plotted in Fig. 7. It is so sharp 
that energy resolution as high as A/ 10^^ has to be used. A nice Lorentz curve appears 
after such an immense blow up. The width of the state in Fig. 6(a) is about 10~^A. 

Fig. 8 is for another set up: the lattice has 14 cells (—10 ~ 3), with a higher barrier 
Vq. It is a cross section of the lowest peak of cell 0. This peak is essentially confined within 
one cell. It is almost as sharp as the peak in Fig. 7(b) even though we are considering a 
relatively small system — there are only 4 cells on the positive x-axis, because the barriers 
are much higher. Its width will be calculated by wave tunneling in Sec. 4. 

The level width will decrease as N increases, because it is more difficult for these 
resonant states to tunnel out. An interesting question is: when our system becomes 
macroscopically large, will the width shrink to zero or to a finite value? This is answered 
in the next section. 

IV. Tunneling and level width 
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A. General formula 

We will be interested in those resonant levels mainly located in a unit cell, say the 
lowest level in cell 0. In the limit of small level width, we may write the level width as 

5E = {n/2) I To p IV (7). 

where v is the attempting frequency in cell 0, and Tq is the transmission coefficient from cell 
to the right edge of the system. It is emphasized that the transmission coefficient used 
here and hereafter refers to the probability currents instead of the probability amplitudes. 
The transmission coefficient through the jth barrier is the same for both directions, and 
is given by 

^ i{kj-i + kj) - K ^ ^ 

On the other hand, the reflection coefficient of a barrier depends on the direction of inci- 
dence, and is given by 

= ^(fcj - kj-i) + K ^-2ikjaj 
^ i{kj + kj-i) - K 

where h refers to reflection from and to the right, and H from and to the left. 

A general formula for Tq may be obtained through a diagrammatic analysis. We 
denote Tj to be the transmission coefficient from the jth cell to the right edge, with the 
understanding that all the paths are contained between the jth cell and the right edge. It 
is then easy to show that 

To = tiTi + tiEiTi + tiE^Ti + . . . 

ti ^ (10) 



1-Ei 

where Ei is the sum over contributions of 'self-energy' loops each of which goes from cell 1 
to the right and resumes the initial position and direction of propagation only in the flnal 
step. See Fig. 9 (a). 

The above formula can be iterated, yielding 



where T,j is the total contribution of 'self-energy' loops starting and ending in cell j. It is 
understood that all the paths for T,j are on the right of cell j, except for the first and the 
last steps which are in the jth cell. 

The formula (11) is an exact result, but the 'self-energy' contribution cannot be ob- 
tained exactly. In the case of strong barriers, we may approximate Ej by rj\-rj+i-\, coming 
from a path reflected once by the {j -|- l)th barrier and once by the jth. barrier. All the 
other paths must have at least two transmissions through the {j + l)th barrier and are 
therefore of higher order. We have done a calculation for the case of A^^, = 4, A = 0.2, and 
Vo = 8. The half width is foimd to be 6.64 x 10~^° by the above approximation, while the 



exact numerical result as shown in Fig. 8 is (6.4 ± 0.1) x 10 . The error is less than 



/o. 



When the barriers to the right of the jth. cell are not strong, we cannot use the above 
approximation. However, if the reflection coeflBcients are small, we only need to sum over 
those loops containing 2 reflections ( Fig. 9 (b) ). Therefore, 



.(2) 



Sj- ~ = J2 tjin^tijTj^ (12) 

where tji = tj^\ • • ■ for j < z, tij = ti - • ■ tj, tjN = tj+i • ■ • tjv- 

T^T'=Y[-^ (13) 

It is interesting to note that the above formula is even valid in the opposite limit of small fs, 
because then only the i = j + 1 term survives. This result is identical to the approximate 
formula that we derived before for the strong barrier case. 

B. 5— function barriers. 

We are now ready to calculate the energy level width of a 5 potential chain subjected 
to electric fleld. When kj ^ k, which is true for a large j such that jA + E :$>Vq, 

~ (14) 

./iTWo/WoTW+m 

\rj\-\ ~ r,-H ~ — , 15) 



If we are considering the limit A — > 0, (14), (15) remain valid if E ^ Vq. 

Using these approximations, we will find that : (phase factors are included) 



(2) 



a 



E 



(16) 



where ^ = (47r/3)-\/A, /? = 1/2 + a/A, a = (7rVo/2)^. This expression is correct as long 
as NA is larger than, or of the same order of magnitude as E. 

The problem now is how to sum up this series? If A is not too small, the difi'erence 

5(^/^/^) ~ (3/2)^/^/^ = 2nVlA ^ 2n. In this case ^/^/^ mod 2n can be considered random. 

(2) 

Therefore, ( is simply written as aj, and Nr is written as N in the following equations 

) ^ ^ 

' ^ (17) 



where ( ) means taking the average with respect to the angle of the complex number. We 
assumed N is so large that the term {N + ^)~^/^ is negligible. 

The quantity \aj\ is small since we are considering j > J. The infinite product 

n(l — aj) ^ in (13) can then be approximated by e ^ when the index j is larger than 
J. The phase of aj is , again, a random number. Therefore, we get 



AT 



N 
3>J 



a / ,iVA + E, 
2A JAT^) 



1/2 



(18) 



with a random phase factor e*^. 
Finally, we get 

N N 



\T'\ = Cj 



j>j j>j 

I N 

Cj exp In \tj I 



exp 



N 



a. 



(19) 



.a NA + E a NA + E 

~ Cjexp I -— ln(— — — ) + W — ln(— — — )cos(^^nd) 



2A ' JA + E 



2 A ' JA + E 



in which Cj is a finite number. We can see that T' will approach zero as AT — > oo 



T' is a good approximation to correct T, at least up to the first order of approx- 
imation. It is improbable that higher order corrections will lead to a nonzero T since 
that would require the correction of the part of ]^ (1 — cTj)"^ diverges at least as fast as 
exp((Q;/2A)ln(A'"/J)). Notice that the dominant term — (q;/2A) In (^ "jA+g ) comes solely 
from the zeroth order contribution. The wavefunction '^{x) ~ is power law local- 

ized, which is a marginal case between extended and exponentially localized cases. Similar 
studies of localization in the energy space also showed marginal behavior [11]. 

C. The Kronig-Penny type barriers 

The same principle can be applied to the Kronig-Penny type barriers. Assume that 
the width and height of the barriers are ai and Vq respectively, that the distance between 
barriers is a2, that a = oi -|- 02 is the size of unit cell, and that the potential drops one 
A over distance a. Then, a straightforward calculation shows that, when j is much larger 
than one 



(1 + k'Jkjf - (1 - k'Jkjfe^^^i''^ 



(20) 



where k'- = ^2m{E + jA)/h; kj = ^/2m{E + jA - Vo)/n. 
When E + jA:$>Vo, 

k'- 1 Vrs 

^l+o^T^ (21) 



kj 2E + jA 

Equation (20) can be further simplified : 
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4^E^+jA) ) (l-c««(2Mi)) (22) 



Therefore, 



ni^.l - I -(t^' £^ UA + E)^ - ^^^(2fe.-^i)) I (23) 

The exponent is a finite quantity. Therefore, up to zeroth order, the level width is finite 
as — > 00. This is true as long as the height of barriers is finite. 

The summation in the exponent of (23) can be approximated by an integral. It turns 
out that, 

nit,|~exp{-a(^y^^^i^| (24) 



4 J A{JA + E) 
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where a is a number of order 1. We can see that, as A becomes very small ( but not smaller 
than E/N ), the level width shrinks to zero with a rate ~ exp(— 1/A). In the A ^ limit, 

(24) is no longer valid. In this case, j A + £" ~ i?, and the tunneling coefficient turns out 
to be ~ exp{—a{Vo/4:)'^{N/E'^)), which is finite as long as the system is finite. 

We can also reproduce the leading term in the exponent of (19) for the case of the 
Kronig-Penny model in the limit that ai — > 0, Vq — > oo, and at the same time keep oiVq 
a constant. 

V. Conclusion 

We have calculated the local density of states p{x,E) for a ID model using delta 
function potential barriers to simulate a periodic potential and a staircase potential to 
simulate an external electric field. The W-S ladders are revealed most clearly as peaks of 
the local density of states in 3D plots of p{x,E). We have demonstrated that the ladder 
structures can be easily washed out in the total density of states. 

We have also offered a new approach, based on the calculation of leakage current in a 
graded array of potential barriers, to evaluate the level width. [12] This complements the 
usual method based on the picture of Zener tunneling between Bloch bands. [13] It is shown 
that, when the system size N is very large, the level width 5E ~ N^~'^/^^ in the case of 5 
barriers. We choose 5 potential barriers to simplify the derivation of the transfer matrices 
as well as transmission and reflection coefficients. The same approach can be applied to 
nonsingular potentials. We found that, when the potential barriers are rectangular, the 
reflection coeflicient vanishes like 1/x for large x, and that the energy level width will not 
shrink to zero as the size of the system becomes inflnite. When the potential is smoother, 
the reflection coefficient vanishes even faster. The electron wave function can be localized 
only in the case of a 5 potential chain, but even then the localization is found to be only 
marginal; that is, it obeys a power law. 

Some of the general properties of W-S ladders are summarized below, and these should 
also be valid in a real solid. Firstly, the number of sequencies presented in a flnite system 
can be estimated from a tilted-band diagram. It can be controlled either by changing the 
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size of the system, or, more conveniently, by changing the electric field. Secondly, among 
these sequences, there is one that is most pronounced. This is best illustrated in Fig. 3 (d). 
The states of the sequence that are closest to the staircase are sharper and more localized 
spatially. This can be explained by the theory in Sec. 4. The states in this sequence have 
initial barriers which are difficult to tunnel through. The states in the other sequences are 
more remote from the origin, and therefore have higher energy with respect to the ground. 
This makes them easier to tunnel through. Thirdly, within the same sequence, the states 
in cell -1 are narrower than those in cell 0, and those in cell are narrower than those in 
cell l...etc. This is because the state in cell -1 has to go through more barriers in order to 
tunnel out. A similar kind of asymmetry leads to the asymmetrical absorption spectra in 
the experiment that was done by AguUo-Rueda et al [9] . 
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Figure Captions 



Fig.l. The potential of the system described in equations (1),(2). It corresponds to 
the case Nr = Ni = 5. 

Fig.2. The plot of p{E) when electric field is zero, (Vq, A,iV^) = (0.2,0,40). Ni in 
this and in all of the following figures are fixed to be 10. 

Fig. 3 (a),(b). Local density of states p{x,E) plotted as a function of x and E. 
{Vo,A,Nr) = (0.2,0.2,40) 12 cells are plotted (-1~10). Range of energy covered is 3 
A, spatial resolution is a/5, energy resolution is A/ 100. The curves will be smoother with 
higher resolution. 

Fig.4 (a),(b),(c),(d). A series of plots of p{E) for a system with (Vq, A) = (0.2, 0.2), 
Nj. = 40, 60, 100, 160 respectively, (a) is for the system plotted in Fig. 3. Also note the 
difference between (a) and Fig. 2. 

Fig.5 (a),(b). Two plots of p{E) for a system with (Vq, A) = (0.2, 1), Nr = 40, 160 
respectively. The electric field of the systems in this figure is 5 times larger than those in 
Fig. 4. 

Fig.6. Plots of the states at Ei = 0.339063 A, Fig. (a) ; and E2 = 0.6757468519 A, 
Fig.(b), in a system with (V, A,7V^) = (1,0.2,40). 

Fig. 7. A cross section of the peak in cell 3 (x = 3.5a) in Fig. 6(b). The range of 
energy covered is extremely small. Ec = 0.67574685190A, half width 6E = 6 x 10~^^A. 

Fig. 8. A cross section of the lowest peak in cell 0. (Vo,A,A^^) = (8,0.2,4). Ec = 
4.7607583641A, 5E = 6.4 x lO^^^A. Notice that the strength of the 5 function barriers 
is 8 time large. The number of cells on the positive x axis is only 4. The level lies only 
slightly below the level of an infinitely deep square well, which is 5 A. 

Fig. 9. (a) A path that makes its last return at point j. We call the bold-faced part a 
self energy loop of j. T,j is the summation of every possible self energy loops of j. (b) A 
self energy loop with only two reflections. The associated tunneling amplitude is stated in 
(12). cr^-^^ is the summation of every possible self energy loops of this kind. 
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